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O ■ Abstract 



Here we generalized a previous model of gravastar consisted of an internal de Sitter spacetime, 



a dynamical infinitely thin shell with an equation of state, but now we consider an external de 

O ' 

Sitter-Schwarzschild spacetime. We have shown explicitly that the final output can be a black hole, 
a "bounded excursion" stable gravastar, a stable gravastar, or a de Sitter spacetime, depending 
on the total mass of the system, the cosmological constants, the equation of state of the thin 
shell and the initial position of the dynamical shell. We have found that the exterior cosmological 
constant imposes a limit to the gravastar formation, i.e., the exterior cosmological constant must be 
smaller than the interior cosmological constant. Besides, we have also shown that, in the particular 
case where the Schwarzschild mass vanishes, no stable gravastar can be formed, but we still have 
formation of black hole. 
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I. INTRODUCTION 

As alternatives to black holes, gravastars have received some attention recently 
partially due to the tight connection between the cosmological constant and a currently 
accelerating universe 3J, although very strict observational constraints on the existence of 
such stars may exist 4]. 

n 

The pioneer model of gravastar was proposed by Mazur and Mottola (MM) [5|. After 
this work, Visser and Wiltshire (VW) [6| pointed out that there are two different types of 
stable gravastars which are stable gravastars and "bounded excursion" gravastars. The first 
one represents a stable structure already formed, while the second one is a system with a 
shell which oscillates around a equilibrium position which can loose energy and to stabilize 
at the end. 

Recently we have done an extensive study on the problem of the stability of gravastars. 
The first model [7| consisted of an internal de Sitter spacetime, a dynamical infinitely thin 
shell of stiff fluid, and an external Schwarzschild spacetime, as proposed by VW [6]. We 
have shown explicitly that the final output can be a black hole, a "bounded excursion" 
stable gravastar, a Minkowski, or a de Sitter spacetime, depending on the total mass m of 
the system, the cosmological constant A, and the initial position Rq of the dynamical shell. 
Therefore, we have shown, for the first time in the literature, that although it does exist 
a region of the space of the initial parameters where it is always formed stable gravastars, 
it still exists a large region of this space where we can find black hole formation. Then, 
we conclude that gravastar is not an alternative model to black hole as it was originally 

n 

proposed by VW models [6j. 

In the second paper jjj, we have generalized the previous work on the problem of stable 
gravastars considering an equation of state p = (1 — 7)0" for the shell, instead of only using 
a stiff fluid (7 = 0). We have found that stable gravastars can be formed even for 7^0, 
since 7 < 1, generalizing the gravastar models proposed until now. We also have confirmed 
the previous results, i.e., that both gravastars and black holes can be formed, depending on 
the initial parameters. 

In the third work we have generalized the former one considering now an interior 
constituted by an anisotropic dark energy fluid. We have again confirmed the previous 
results, i.e., that both gravastars and black holes can be formed, depending on the initial 
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parameters. It is remarkable that for this case we have an interior fulfilled by a physical 
matter, instead of a de Sitter vacuum. Thus, it is similar to phantom energy star models. 

Nowadays, several kinds of observational data indicate that our universe is in accelerated 
expansion. In Einstein's general relativity, in order to have such an acceleration, one needs 
to introduce a component to the matter distribution of the universe with a large negative 
pressure. This component is usually referred as dark energy. Astronomical observations 
indicate that our universe is flat and currently consists of approximately 2/3 dark energy 
and 1/3 dark matter. The nature of dark energy as well as dark matter is unknown, and 
many radically different models have been proposed, such as, a tiny positive cosmological 
constant. Based on this fact, we would like to ask how the picture of the evolution of 
gravastar formation is influenced by an exterior spacetime with a positive cosmological 
constant. 

Recently, Carter [h]] studied spherically symmetric gravastar solutions which possess an 
(anti) de Sitter interior and a (anti) de Sitter-Schwarzschild or Reissner-Nordstrom exterior. 
He followed the same approach that Visser and Wiltshire took in their work 6] assuming a 
potential V(a) and then founding the equation of state of the shell. He found a wide range 
of parameters which allows stable gravastar solutions, and presented the different qualitative 
behaviors of the equation of state for these parameters. 

Differently from Carter's work [h]], we consider here another approach. We generalize 
our second work in gravastars {s| , introducing an external de Sitter-Schwarzschild spacetime, 
to study how the cosmological constant affects the gravastar formation. We first assumed an 
equation of state, p — (1 — 7)p, and, using Israel conditions, derived a potential depending 
on the parameters of the interior, the shell and the exterior of the gravastar's prototype. 
We, then, studied the types of compact objects that can be generated according to this 
potential, to the parameters related to the cosmological constants and to the masses of our 
model. We found that both gravastars and black holes can be formed. 

The paper is organized as follows: In Sec. II we present the metrics of the interior and 
exterior spacetimes, with theirs extrinsic curvatures, the equation of motion of the shell 
and the potential of the system. In Sec. Ill we discuss the particular cases where the 
Schwarzschild mass is null, and another where we have the same cosmological constant in 
the interior and the exterior of the thin shell, which is presented in section IV. In Sec. V 
we investigate the formation of gravastar from numerical analysis of the general potential. 
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Finally, in Sec. VI we present our conclusions. 



II. FORMATION OF GRAVASTARS IN A DE SITTER-SCHWARZSCHILD 
SPACETIME 



The interior spacetime is described by the de Sitter metric given by 



dsj = -hdt 2 + f 2 dr 2 + r 2 dtt 2 , 



(1) 



where f\ = 1 - (r/L t ) 2 , f 2 = 1 _ (r ? /£ . )3 , U = ^3/A t and dQ 2 = d6 2 + sin 2 (fl)d0 2 
The exterior spacetime is given by a de Sitter-Schwarzschild metric 



ds z e = -fdv 2 + f^dr 2 + rW, 



where f = \- 2 -f- (r/L e ) 2 and L e = ^3/A~ e . 

The metric of the hypersurface do the shell is given by 



(2) 



ds% = -dr 2 + R 2 (r)dQ 2 



(3) 



Since ds 2 = ds 2 = ds| then r s = r s = R, and besides 

-l 



and 



t 



-a 



/.-/.if 



5 2m /i? 

\T e . 



R 2 



where the dot represents the differentiation with respect to r. 
Thus, the interior and exterior normal vector are given by 



2m ( R 
~R ~ 



(4) 



(5) 



< = (-i?,t,0,0), 



(6) 
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and 



n;. 



(-R,v,0,0). 

The interior and exterior extrinsic curvature are given by 



(?) 



Kl 



[(3LfR - L\t 2 + 2LjR 2 t 2 - Rt 2 )Rt - (U + R)(Li - R)(Rt - Rt)Lf] x 



{L l + R)-\L l -R)- l L^ 



Kl e = i(U + R)(U-R)LfR 



(8) 

(9) 
(10) 



K e TT = v\(2L\mv + L\RR - L 2 Rv + R 3 v)(2L 2 e mv - L 2 RR - L 2 e Rv + R 3 v) - 

2L A e R 2 R 2 }((2m - R)L\ + R 3 )-\L 2 e m - R 3 )L^R~ 3 + Rd - Rv (11) 



3^ r -2 

e 



K% e = -v((2m - R)L\ + R 3 )L 
K e u = K e ee sin 2 (9). 



Since we have 
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[K ee ] = K e ee - Kl e = —M, 
where M is the mass of the shell, thus 
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M = vR 



tR 



1- ^ 



Then, substituting equations (TjJ and ([5]) into ( |T5l) we get 



M + R 



1/2 



— i? 



1-1^) + # 



1/2 



0. 



(12) 
(13) 

(14) 
(15) 



(16) 



Solving the equation fflBT) for R 2 /2 we obtain the potential V(R,m, Li, L e ). In order to 
keep the ideas of our work [8j as much as possible, we consider the thin shell as consisting 
of a fluid with a equation of state, a = (1 — 7)1?, where a and $ denote, respectively, the 
surface energy density and pressure of the shell and 7 is a constant. The equation of motion 



of the shell is given by 



111 



M + teRRd = AnR'lTap^n 13 } = irR 2 [T+u%n1 - T~m 



? 2 [rp+ a^P 



a 13 

■ n_ 



(17) 
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where u a is the four-velocity. Since the interior and the exterior spacetimes correspond to 
vacuum solutions, we get 

M + 8irRR(l - 7)0- = 0, (18) 
and since a = Mj (AtcR 2 ) we can solve equation fjlSp giving 



M = kR 2 ^. 



(19) 



where k is an integration constant. 

Substituting equation f[T9"j) into V(R,m, L i: L e ) we obtain 

V(R,m,Li,L e ,k,j) = 
1 



8R 2 UUk 2 



[-AR 2 L 4 e Ltk 2 + ARmLtLtk 2 + 2R^L\k 2 L\ + R^^L\ 



-4R(- 4 ^LjL 2 m - 2R^ +12 ^L 2 e L 2 + 2R 4 L 4 e L 2 k 2 + h\m 2 h\ 

+4i?( _4 7+ 9) L 4 mL 2 + #(-47+12) £4 + #(47-4) £4^4] . 



(20) 



Redefining the Schwarzschild mass m, the cosmological constants Lj and L e and the 
radius R as 



m = mk 2 t- 3 . 



Lik^-z. 



L P = L P k 2 i-z, 



R = Rk~^. 



we get the potential 



V(R,m,L i} L e ,j) = 

+ R + 4 + K + 2L 2 

#(-47+10) #2 m #(-47+7) #(-47+10)- 
+ + 



2L?L 2 



2L 2 



L 2 



4L! 



Redefining L e = aLj we finally get 



m #(-47+7) #(-4 7 +10) 
+ 



4Lf 



(21) 
(22) 
(23) 
(24) 



(25) 



V{R, m, a, L h 7) = 

i? 4 2L 2 



m #(-47+7) #(-4 7 +10) 
+ 



4L? 



#(-4 7 +io) 



+ 



i? 2 



2a 2 L? 2a 2 L 2 



+ 



#(-47+7) #(-47+10) 



a 



2 L 2 



+ 



4a 4 Lf j 



(26) 
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TABLE I: This table summarizes the matter classification based on the energy conditions of the 
shell, in terms of the parameter 7. 



Matter 



Condition 1 



Condition 2 



7 



Standard Energy 

Dark Energy 
Phantom Energy 



a + 2p > 
a + 2p < 
a + 2p < 



a + p > 
a + p > 
a + p < 



7 < 1-5 
1.5 < 7 < 2 

7 > 2 



It is curious to note that this potential is independent of the sign of the parameter a. 

Therefore, for any given constants m, a, Li and 7, equations ( }25l) or (1261) uniquely 
determines the collapse of the prototype gravastar. Depending on the initial value Rq, the 
collapse can form either a black hole, or gravastar, or a de Sitter spacetime. In the last case, 
the thin shell first collapses to a finite non-zero minimal radius and then expands to infinity. 
To guarantee that initially the spacetime does not have any kind of horizons, cosmological 
or event, we must restrict Ro to the ranges simultaneously, 

2m<R < L u (27) 

2m< R < L e , (28) 

where Ro is the initial collapse radius. 

In order to fulfill the energy condition a + 2p > of the shell and assuming that p = 
(1 — 7)er we must have 7 < 1.5. On the other hand, in order to satisfy the condition 
a + p > 0, we get that 7 < 2. The dominant energy condition is only satisfied for < 7 < 2. 
Although the phantom energy is usually considered as a kind of dark energy, in this paper 
we will use the expression dark energy for the case where the condition o +p > is satisfied 
and phantom energy otherwise. Hereinafter, we will use only some particular values of the 
parameter 7 which are analyzed in this work. See Table I. 

Since the potential, equations ( 1251) or ( 1261) . is very complex to manipulate analytically, 
we have analyzed several special cases. 



III. CASE m = 



This case represents a system where the Schwarzschild mass vanishes and the combination 
of both cosmological constant (interior and exterior) imposes a very special junction thin 
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shell. Note that from equation (fT6|) . this configuration is possible only if a ^ 1, otherwise if 
a = 1 then we have M = 0, i.e., the thin shell vanishes. 
From the equation (|2"oT) we get 



1 #( 4 T~ 

2 8~ 



R 2 R (-i 1+ lQ) ^(_4 7 +10) 



i? 2 



4L 2 8L 4 4a 2 L 4 4a 2 L 2 

and differentiating the potential dV(R)/dR, we get 



^(-47+10) 

8a 4 L 4 







(29) 



#(47-6)(4 7 _ 6 ) ^ jR (-4 7 +10)(_4 7 + l0) ^ _R(-47+10)(_4 7+10 ) 



8i? 2L 2 
R J R(- 4 ^ +1 °)(-47 + 10) 



8i?L 4 



4Ra 2 L 4 



0. 



(30) 



2a 2 L 2 8i?a 4 L 4 
From these two equations we can obtain the point where the potential has a minimum and 
equal to zero. Solving simultaneously the equations fl29|) and fl30|) we get 

2 



2"2^3 



{-2a 4 + a 4 7 - 2a 2 7 + 6a 2 + 7 - 2+ 



a 2 (-2 + 7 ) 

[(a 2 + l) 2 (4a 4 + a 4 7 2 - 4a 4 7 - 7a 2 - 2a V + 8a 2 7 + 4 + 7 2 - 4 7 )] * } } 



2(27-3) 



(31) 



{-(-4 + R A ~<~ & ) [R 2 a 2 + R 2 + (R 4 a 4 + 2R 4 a 2 + R 4 + AR^ +10 ^ 



(-4 + R^~ Q )a 

gjR (-47+10) a 2 + 4jR (-47+10) a 4 _ ^(4 7 -6) £,(-47+10) + 2 R^~^ R^ +1 ^ 0? - R^+W a 4 )' 



For a = oo we get the same results of previous work 

2(7-2 



8j, given by 



R r 



2 7 -5 



i 

3-2 7 



and 



2 7 - 5 



1/2 



(33) 



(34) 



2 7 -3 

We can see from figure CD that the quantity V"(R), calculated at the values R = R c and 
Li = L l c , is always negative, for a large range of values for a and 7 ( — 1 < 7 < 1.5 and 
< a < 100). This means that, if we impose V(R) = V'(R) = 0, we have always formation 
of black holes, instead of formation of stable gravastars. 

In the next sections, we will analyze another interesting particular case, where m 7^ 
and a — 1. 



(32) 
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0.5 
gamma 



FIG. 1: This plot shows, in terms of a and 7, the second derivative of the potential V(R, m, a, Li, 7) 
with respect to R, calculated at the values R = R c and Lj = L*, in the intervals —1.5 < 7 < 1.5 
and < a < 100, for m = 0. 

IV. CASE m + AND a = 1 

In this case we consider L e = Li = L since a = 1. 
From the equation we get 

-1 + ^ + 7 i? 4y ~ 6 + m 2 iT 4 ^ 4 + = (35) 
R 4 L 

and differentiating the potential dV(R)/dR, we get 

m 1 i? 47 - 6 (47 -6) m 2 iT 47+4 (-4 7 + 4) R , , 
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From these two equations we can obtain the point where the potential has a minimum and 
equal to zero. Solving simultaneously the equations fl35|) and fl36|) we get 

(-3 + v /(25 + 16 7 2 i? 4 ^ - 40 ^R^ + 32 Wy - 16 i? 6 ) i?" 4 ^) 4i? 47 

m < = " IFW^T) " • < 37 > 



L 



c 



-8 J R 13 - 8 ^(2 7 - if 



16 jR n-8 7 (3 7 _ 2 7 2 - 1) + (47 - 5)i?- 4 ^+ 5 [47 - 5 + ^(47 - 5) 2 + 16 #-^+6(27 - 1)] ' 

(38) 

From the figures [2] and [3] we can see that the critical mass m c is positive only in the range 
0.5 < 7 < 1.5. Besides, from the figure H]we can note that there is not any real value for the 
critical cosmo logical constant L c in the interval 7 < 1. As a consequence of these results, 
the second derivative of the potential V" (R), shown in the figure El is negative for 7 < 1 and 
positive for 1 < 7 < 1.5. For 7 = 1 we have V" (R) = implying that we have an inflection 
point in the potential. 

Combining all these facts, we conclude that for L e = Li we obtain the following: 

1. For 7 > 1.5, which corresponds to a dark energy shell, none structure is formed. 

2. For 0.5 < 7 < 1.5, which corresponds to a standard fluid shell, it can collapse to 
a black hole (0.5 < 7 < 1), or it does not collapse, reaching an equilibrium stage, 
forming a stable gravastar (1 < 7 < 1.5). 

3. For 7 = 1 none gravastar is formed. 

Then, for a = 1, we have shown that no stable gravastar can exist, for 7 < 0.5. 

V. GENERAL CASE 

The expressions for the potentials in the present case makes difficult a complete analytic 
analysis, so we shall study it numerically. Our main strategy is to start with the values of 
m c obtained for the case studied in our previous work 8J], where m ^ and Li ^ 00, and 
then gradually turn on L e . The potential V(R, m c , L i: L e , 7) is plotted as a function of R, 
by finely tuning L e until a stable gravastar or a "bounded excursion" gravastar is found. We 

11 



FIG. 2: This plot shows that, in terms of R and 7 and for L e = Li, the critical mass m c is always 
negative, when V(R) = and V'(R) = 0, for the interval 0.5 < 7 < 1.5, since the critical mass is 
not defined for 7 < 0.5. 

also made another approach, solving the system of equations V(R,m c , Li, L^j) = and 
V'(R,m c , Li, L e ,j) = for R and L e and fixing the parameters m c , Lj and 7 in order to 
compare the results we obtained for L e . It was seen that there is a range of L e in which 
"bounded excursion" stable gravastars are found, i.e., L™ m < L e < L™ ax . For L e > L™ ax 
we have found only stable gravastars. 

We must call attention to the fact that, hereinafter, we will not consider the physical 
situation where there is dispersion of the star. If the initial radius of the collapse is greater 
enough, the star will first contract to its minimal radius and then expand to infinity, whereby 
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FIG. 3: This plot shows that, in terms of R and 7 and for L e = Li, the critical mass m c is always 
negative, when V(R) = and V'(R) = 0, for the interval 1.5 < 7 < 2. 



a de Sitter spacetime is finally formed. 

Figures IBTfTTI show the behavior of the potential as a 



7 = — 1 . This case was not studied on our previous work ] where there was no cosmo logical 
constant external to the thin shell. So, we used the analytic expression (2.23) from our 
previous work Q| to calculate m c . This situation is analogous to our present work if we use 
the potential with Lj = L e = 00. The potential is shown in figure [6] where m c ~ 0.5055981490 
and R c m 1.023836256. For m > m c the potential V(R) is strictly negative and the collapse 
always forms black holes. For m = m c , there are two different possibilities, depending on 
the choice of the initial radius Rq. In particular, if the star begins to collapse with Rq > R c , 



unction of R for the case where 
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200 

r 300 

400 

500 

FIG. 4: This plot shows, in terms of R and 7, the cosmological constants L e = Li when V(R) = 
and V'(R) = 0, for the interval 0.5 < 7 < 1.5, since the cosmological constant is not defined for 
7 < 0.5. The critical mass m c is always negative for 1.5 < 7 < 2.0. 

the collapse will asymptotically approach the minimal radius R c . Once it collapses to this 
point, the shell will stop collapsing and remains there for ever. However, in this case this 
point is unstable and any small perturbations will lead the star either to expand for ever 
and leave behind a flat spacetime, or to collapse until R = 0, whereby a Schwarzschild black 
hole is finally formed. On the other hand, if the star begins to collapse with 2m c < Rq < R c , 
the star will collapse until a black hole is formed. For m < m c , the potential V(R) have a 
positive maximum, and the equation V(R, m < m c ) = has two positive roots R\£ with 
R2 > R\ > 0. There are two possibilities here, depending on the choice of the initial radius 
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FIG. 5: This plot shows, in terms of a and 7, the second derivative of the potential V(R, m,a = 
1, Li, 7) with respect to R, calculated at the values R = R c and Lj = U c , for the interval 0.4 < 7 < 
1.5, since the derivative is not defined for 7 < 0.5. 

Rq. If Rq > R2, the star will first collapse to its minimal radius R = R2 and the expand to 
infinity, whereby a Minkowski spacetime is finally formed. If 2m c < Rq < R\, the star will 
collapse continously until R=0, and a black hole will be finally formed. As we always have 
V" < 0, it means that no stable stars exist in this case. 

For the case of the figure El i.e, m = 0.5055981490, Lj = L e = 00 and 7 = —1, we have 
analyzed the behavior of the potential for the parameter 7 = 0.7, 1.7, 3 and we have found 
that we get only dispersion of the shell. 

The figures [7| and M show the case where L e = 00, but L, 7^ 00. Variations of m fixing the 
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parameter Li and variations of Li fixing the parameter m reveal that both stable gravastars 
and "bounded excursion" stable gravastars can be formed, but not excluding the existence 
of black holes. 

For the general case, where both Lj and L e are not infinity, it is shown the potential 
V(R) as a function of R for some specific values of 7, which are 7 = —1, 7 = 0, 7 = 0.4 
and 7 = 0.7 representing standard energy, 7 = 1.7 representing dark energy and 7 = 3 for 
phantom energy. Note that 7 = — 1 and 7 = 3 violate the dominant energy condition. Note 
also that in the Carter's work [lOj, the dominant energy condition is considered to restrict 
acceptable solutions. In our case this corresponds to the cases 7 = 0, 0.4, 0.7 and 7 = 1.7. 
We found that the shell must have standard energy (figures [HI CE21 EH and [H} in order to 
have both stable gravastars or "bounded excursion" stable gravastars (the later existing 
whenever L™ m < L e < L™ ax as explained in the text), but never excluding the existence of 
black holes or the formation of a de Sitter space depending on the choice of initial radius 
Ro (It is important to verify the restriction on the values Ro can assume, obeying both 
2m < Rq < Li and 2m < Rq < L e .). For dark energy shells and for phantom energy shells 
there are not formation of gravastars (figures fT5l and [TBI) . Variations of 7 in the potentials 
studied also show that when the region of 7 represents dark or phantom energy, there are 
only possibilities of formation of black holes or de Sitter spacetime (figures ED, HH [20] and 
[231 . When the shell is made of standard energy we can have gravastars or black holes (figures 
ITOl [T71 [T9l [2T1 and 1221) . Thus, we only find gravastars for standard energy shells, satisfying 
or not the dominant energy conditions. 



VI. CONCLUSIONS 



In this paper, we have generalized the problem of the stability of gravastars studied 
recently by us [8], introducing a positive cosmological constant in the exterior spacetime. 
Thus, the model consists of a de Sitter interior spacetime, a dynamical infinitely thin shell 
of fluid with an equation of state p = (1 — 7)0", and an external de Sitter-Schwarzschild 
spacetime. We have shown explicitly that the final output can be a black hole, a "bounded 
excursion" stable gravastar, a stable gravastar, or a de Sitter spacetime, depending on the 
total mass m of the system, the parameter a, the constant Li, the parameter of the shell 
7 and the initial position Rq of the dynamical shell. All these possibilities have non-zero 
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TABLE II: This table summarizes all possible kind of energy of the interior fluid and of the shell. M, 
S, dS and dSS denote Minkowski, Schwarzschild, de Sitter and de Sitter-Schwarzschild spacetimes, 
respectively. 



Case 


Interior 


Shell Energy 


Exterior 




Figures 




Conditions 


Structures 


A 


M 


Standard 


S 




6 




m > m c 


Black Hole 


B 


dS 


Standard 


S 




7 




m = m c , m = 0.51 


Gravastar 












7 




m = 0.53 


Black Hole 












8 




U = L c , Li = 1.4 


Gravastar 












8 




L % = 1.0 


Black Hole 








dSS 


9, 


12, 13, 


14 


mini ^ <^ max 


Gravastar 










9, 


12, 13, 


14 


L e < L e min 


Black Hole 












10, 11 




7 = -l 


Gravastar 
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7 = 0.0 


Gravastar 












19, 21 




7 = 0.4 


Gravastar 












19, 21 




7 = 0.0 


Black Hole 












22 




7 = 0.7 


Gravastar 












22 




7 = 0.0,0.4 


Black Hole 


C 


dS 


Dark 


dSS 
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m > m c 


Black Hole 


D 


dS 


Phantom 


dSS 


11. 


, 18, 20, 


23 


7 = 3 


Black Hole 
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m > ra c 


Black Hole 



measurements in the parameter space of m, Li, a, 7 and Rq, for both gravastar and black 
hole. 

For m — 0, the analysis of the potential has shown that, if we impose V(R) = V'(R) = 0, 
we have always formation of black holes, instead of formation of stable gravastars. Compar- 
ing the results from [8] (L e = 00) with this work, we have confirmed that, in a more general 
way, there is no formation of gravastar even with the introduction of a L e ^ 00. 

On the other hand, for m 7^ 0, if Lj = L e (a = 1) we have formation of black hole 
or stable gravastar. These gravastars are only possible for 1 < 7 < 1.5, satisfying all the 
dominant energy conditions. It is interesting to remark that this case can not be compared 
to other one already studied by us [§]], except for Li = L e = 00, which was shown in the 
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figures 6 and 7, in that paper, and in the figure 6 of this work. While we have gravastars 
there for 7 < 1, here the gravastar formation is limited to 1 < 7 < 1.5, showing that these 
intervals are complementary to each other, except for 7 = 1. 

In the general case, i.e., m ^ 0, Lj ^ L e , it was seen that there is a range of L e in which 
"bounded excursion" stable gravastars are found, i.e., L™ m < L e < L™ ax . (Reminding that 
the curve for L e = 00 is very close to the curve for L e = L™ ax .) Stable gravastars were 
found for L e > ]j^ ax . Besides, this interval depends on the values of Lj_and 7. Let us now 
compare figures 12 and 13 , presented here, with figures 8 and 10, from 8J, respectively. We 
can state that, from figures 8 and 10 [S], the bigger is Li (for L e = 00) the bigger is the 
tendency to the collapse of the shell, forming a "bounded excursion" gravastar or a black 
hole. Moreover, from figures 9, 12, 13 and 14 of this paper, for a given L i: the formation of 
gravastars depends on the value of L e (L e > L e mm , with L e mm > Li) in a such way that, 
instead of what occurs for Lj, the smaller is L e the bigger is the tendency to the collapse. 
These conclusions are in agreement to the gravastar requirement proposed by Horvat & Ilijic 
[l]. The reason is that the dark energy density inside the gravastar have to be greater than 
the surround spacetime, i.e., Li < L e . All these results can be summarized in Table II. 
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FIG. 6: The potential V(R) for 7 = -1, m c = 0.5055981490, L e = 00, and U = 00 (the second 
curve top-down). The others curves represent values for m < m c (first curve top-down) and 
m > mc (the third and fourth curve top-down). Case A 
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FIG. 7: The potential V(R) for 7 = -1, m c = 0.5055981490, U = L C = 1.932841686 and L e = 00 
(second curve top-down). The first curve top-down assumes m = 0.49. The third and fourth curves 
top-down assume m = 0.51 and m = 0.53, respectively. Case B 
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FIG. 11: The potential V(R) for 7 = -1, m c = 0.5055981490, U = 1.932841686 and Lf ax = 
468695.8789. (the curve that has a minimum). The others two curves top-down use the values 
7 = 1.7 and 7 = 3, respectively. Cases B and D 
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FIG. 12: The potential V{R) for 7 = 0, L e = 133843.0443, U = 2.8743397865 and 
m c = 0.5170643255 (the first curve top-down). The second curve top-down is calculated using 
6.16479753 < L e < 133843.0443. The third curve top-down is obtained assuming L e = 6.16479753. 
The fourth curve top-down assumes L e < 6.164749753. The curve for L e = 00 is very close to the 
curve for L e = L™ ax . These curves generalize the results presented in the figure 8 from 8J]. Case 
B 
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FIG. 13: The potential V(R) for 7 = 0.4, L e = 164286.4284, U = 4.3336020545 and 
m c = 0.5350989795 (the first curve top-down). The second curve top-down is calculated us- 
ing 7.828275741 < L e < 164286.4284. The third curve top-down is obtained assuming L e = 
7.828275741. The fourth curve top-down assumes L e < 7.828275741. The curve for L e = 00 is 
very close to the curve for L e = L™ ax . These curves generalize the results presented in the figure 
10 from \8\. Case B 
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FIG. 14: The potential V(R) for 7 = 0.7, L e = 127780.9030, Li = 9.1891232 and m c = 
0.5794043665 (the first curve top-down). The second curve top-down is calculated using 
13.53811906 < L e < 127780.9030. The third curve top-down is obtained assuming L e = 
13.53811906. The fourth curve top-down assumes L e < 13.53811906. The curve for L e = 00 
is very close to the curve for L e = L™ ax . These curves generalize the results presented in the figure 
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FIG. 15: The potential V(R) for 7 = 1.7, U = 0.6064576241 x 10 8 and m c = 1.043246242 (the first 
curve top-down). The others curves represent m > m c . The potential is insensible for variations 
of L e . These curves generalize the results presented in the figure 20 from 8]. Case C 
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FIG. 16: The potential V(R) for 7 = 3, L e = 549993006.6, U = 10410.51705 and m c = 
0.5120894280 (the first curve top-down). The second curve represent m > m c . Case D 
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FIG. 17: The potential V(R) for 7 = 0.7, L e = 133843.0443, U = 2.8743397865 and m c = 
0.5170643255 (the first curve top-down). The second curve top-down is calculated using 7 = 0.4. 
The third curve top-down is obtained assuming 7 = 0. Case B 
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FIG. 18: The potential V(R) for 7 = 0, L e = 133843.0443, U = 2.8743397865 and m c = 
0.5170643255 (the first curve top-down). The second curve top-down is calculated using 7 = 1.7. 
The third curve top-down is obtained assuming 7 = 3. Case D 



31 




FIG. 19: The potential V(R) for 7 = 0.7, L e = 164286.4284, U = 4.3336020545 and m c = 
0.5350989795 (the first curve top-down). The second curve top-down is calculated using 7 = 0.4. 
The third curve top-down is obtained assuming 7 = 0. Case B 
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FIG. 20: The potential V(R) for 7 = 0.4, L e = 164286.4284, U = 4.3336020545 and m c = 
0.5350989795 (the first curve top-down). The second curve top-down is calculated using 7 = 1.7. 
The third curve top-down is obtained assuming 7 = 3. Case D 
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FIG. 21: The potential V(R) for 7 = 0.7, L e = 127780.9030, L; = 9.1891232 and m c = 
0.5794043665 (the first curve top-down). The second curve top-down is calculated using 7 = 0.4. 
The third curve top-down is obtained assuming 7 = 0. Case B 
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FIG. 22: The potential V(R) for 7 = 0.7, L e = 127780.9030, L; = 9.1891232 and m c = 
0.5794043665 (the first curve top-down). The second curve top-down is calculated using 7 = 0.4. 
The third curve top-down is obtained assuming 7 = 0. Case B 
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FIG. 23: The potential V(R) for 7 = 0.7, L e = 127780.9030, L; = 9.1891232 and m c = 
0.5794043665 (the first curve top-down). The second curve top-down is calculated using 7 = 1.7. 
The third curve top-down is obtained assuming 7 = 3. Case D 
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